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Abstract 

The inclusive production of D* ± (2010) mesons in deep-inelastic e ± p scattering is 
measured in the kinematic region of photon virtuality 100 < Q 2 < 1000 GeV 2 and in- 
elasticity 0.02 < y < 0.7. Single and double differential cross sections for inclusive 
D* meson production are measured in the visible range defined by \rj(D*)\ < 1.5 and 
Pt{D*) > 1.5 GeV. The data were collected by the HI experiment during the period from 
2004 to 2007 and correspond to an integrated luminosity of 351 pb -1 . The charm con- 
tribution, F^ , to the proton structure function F2 is determined. The measurements are 
compared with QCD predictions. 
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1 Introduction 



The measurement of the charm quark production cross section in deep inelastic scattering (DIS) 
at HERA is a powerful means of testing perturbative quantum chromodynamics (QCD). Within 
this framework, a significant contribution to charm production arises from the boson-gluon 
fusion process which is sensitive to the gluon density in the proton. With increasing photon 
virtuality, Q 2 , the charm contribution to the inclusive ep scattering cross section rises from a few 
to up to 20%. Therefore, the treatment of the effects related to the charm quark contribution, 
in particular the mass effects, in perturbative QCD calculations is an important issue in the 
determination of parton distribution functions (PDFs). Different schemes to incorporate these 
effects are available. 

Previous measurements were performed by identifying charm quarks via D mesons flU 13 
or using variables which are sensitive to the lifetime of heavy flavour hadrons 0. This 
paper presents a measurement of the D*^ 1 meson production cross section in the range of large 
photon virtualities 100 < Q 2 < 1000 GeV 2 . The data were collected with the HI detector at 
HERA during the running period 2004 — 2007 when HERA operated with 27.6 GeV electrons^ 
and 920 GeV protons colliding at a centre of mass energy of y/s = 319 GeV and correspond 
to the integrated luminosity of 351 pb" 1 . The measured cross sections are compared to QCD 
predictions providing an insight into the dynamics of D* ± meson production at high Q 2 . The 
charm contribution, F 2 CC , to the proton structure function F 2 is determined. 



2 HI Detector 

A detailed description of the HI detector can be found elsewhere Q. In the following only 
detector components relevant to this analysis are discussed. A right handed coordinate system 
is employed with the origin at the position of the nominal interaction point that has its z-axis 
pointing in the proton beam, or forward, direction and x(y) pointing in the horizontal (vertical) 
direction. The pseudorapidity is related to the polar angle 9 by r\ = — lntan(6 l /2). 

Charged particle tracks are reconstructed in the central tracking detector (CTD). It con- 
sists of two cylindrical central jet drift chambers (CJC) placed concentrically around the beam- 
line, complemented by the silicon vertex detector |[6l, inside a solenoid with a homogeneous 
magnetic field of 1.16 T. The CJCs are separated by a drift chamber which improves the in- 
coordinate reconstruction. A multiwire proportional chamber mainly used for triggering Q 
is situated inside the inner CJC. The CTD provides a particle momentum measurement over 
the polar angle 15° < 9 < 165°. The trajectories of charged particles are measured with 
a transverse momentum resolution of ct(pt)/pt ~ 0.002 p T /GeV © 0.015. The interaction 
vertex is reconstructed from CTD tracks. The Liquid Argon (LAr) calorimeter [[8]| is used to 
measure the energy and direction of electrons, photons and hadrons. It covers the polar angle 
range 4° < 9 < 154° with full azimuthal acceptance. Electromagnetic shower energies are 
measured with a precision of a(E)/E = 12%/^/E/GeV © 1% and hadronic energies with 

'in this paper "electron" is used to denote both electron and positron. 



4 



a(E)/E = 50%/ y/E/ GeV © 2%, as determined in test beam measurements [9|. In the back- 
ward region, energy measurements are provided by a lead/scintillating-fibre (SpaCal) calorime- 
ter [[T0]| covering the angular range 155° < 9 < 178°. For electrons a relative energy resolution 
of a(E)/E = 7%/ ^/E/GeV © 1% is reached, as determined in test beam measurements [fTTTl . 
The SpaCal also provides time-of-flight information for trigger purposes. The luminosity is 
determined from the rate of the Bethe-Heitler reaction ep — > epj, measured using a photon 
detector located close to the beam pipe at z = — 103 m, in the backward direction. 



3 Models of Open Charm Production 

Open charm production in electron-proton collisions can be described within different schemes. 
At energy scales larger than the charm quark mass, calculations can be performed within the 
zero-mass variable-flavour-number scheme (ZMVFNS) lfT2ll . where the charm quark is treated 
as a massless parton in the proton. The fixed-flavour-number scheme (FFNS) [13] applies close 
to the charm production threshold and takes into account heavy quark mass effects. In the 
latter scheme all quark flavours lighter than charm are treated as massless with massive charm 
being produced dynamically via boson-gluon fusion. A consistent treatment of heavy quarks in 
perturbative QCD over the full energy scale range should be provided through the generalised 
mass variable flavour number scheme (GMVFNS) lfl~4"l . 

The prediction of open charm production in FFNS at next-to-leading order (NLO) uses 
separate programs to calculate inclusive lfT3l and exclusive [[15] (HVQDIS) quantities. The 
momentum densities of the three light quarks and the gluon in the proton are evolved using 
the DGLAP equations OH. For the proton structure the FFNS PDF set MRST2004FF3 El is 
used. The charm quark mass is fixed to m c = 1.43 GeV in accordance with this PDF set. The 
renormalisation and factorisation scales are set to [i r = /j,f = // = \/Q 2 + Am 2 .. The charm 
fragmentation fraction into D*^ mesons is taken as f(c — > D*) = 23.8 ± 0.8% lfi~8l from the 
combination of measurements in e + e _ experiments. 

In the ZMVFNS calculation at NLO lfi~2ll a charm mass of 1.6 GeV, renormalisation and 
factorisation scales of fi r = fif = /io = ^Q 2 + Am 2 c and the CTEQ6.6M QH parton densities 
are used. The perturbative fragmentation function |[20l is evolved to the chosen scale of the 
transverse D*^ 1 momentum in the photon-proton rest frame, p^(D*). 

Events containing charm quarks are generated using the Monte Carlo programs RAPGAP GTI 
and CASCADE [22] and are passed through a detailed simulation of the detector response to 
determine the acceptance and efficiency and to evaluate the systematic uncertainties associated 
with the measurements. 

The RAPGAP program, based on collinear factorisation and DGLAP evolution, is used to 
generate events containing cc pairs via photon-gluon fusion. The leading order (LO) matrix 
element with massive charm quarks is used. Parton showers, based on the DGLAP evolution, 
model the higher order QCD effects. The charm quark mass is set to 1.43 GeV. The proton 
structure is described by the PDF set CTEQ6.5M [23] and the factorisation and renormalisation 
scales are set to /j, r — /if — [Iq — \JQ 2 + p\. 
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The CASCADE program is based on the k T factorisation approach. This calculation of 
the photon-gluon fusion matrix element takes into account the charm quark mass as well as 
the virtuality and transverse momentum of the incoming gluon. Gluon radiation from the in- 
coming gluon as well as parton showers from the outgoing charm and anti-charm quarks are 
implemented in a manner which includes angular ordering constraints. The gluon density of the 
proton is evolved according to the CCFM equations E4l . The charm quark mass and the renor- 
malisation scale are set to m c = 1.5 GeV and (i r = \JQ 2 + pj>, respectively. The unintegrated 
gluon distribution is described by the parametrisation set AO [|25l . 

The kinematics of production depend not only on the charm quark production but 
also on the c — > -D* ± fragmentation process. The charm fragmentation function has been 
measured at HI E6l using inclusive D*^ 1 meson production. The Kartvelishvili fragmentation 
function 112711 . which is controlled by a single parameter a, is used. The parameter values cor- 
responding to the programs used in the present analysis are shown in Table [Q They depend on 
the centre of mass energy squared of the hard process, s. To obtain the visible D*^ 1 production 
cross sections in HVQDIS, charm quarks are fragmented independently in the photon-proton 
centre of mass frame into D*^ mesons according to Kartvelishvili function. In the RAPGAP 
and CASCADE programs hadronisation is performed using the Lund String Model ll28l |29ll . 
The momentum fraction of the charm quark carried by the D*^ 1 meson is modelled according 
to the Bowler parameterisation [|30l . The longitudinal part of the fragmentation function is 
reweighted to the Kartvelishvili function. 



Model 


s < 70 GeV 2 


§ > 70 GeV 2 


HVQDIS 
RAPGAP 
CASCADE 


a= 6.0±\i 
a = 10.3lJ;g 
a= 8At\i 


a = 3.3j£t 
a = 4.4±8;| 
« = 4.5^ 



Table 1: Parameter a of the Kartvelishvili fragmentation function as used in the analysis. 



The contribution of beauty production is estimated using the HVQDIS calculation, with 
hadronisation corrections determined using RAPGAP. The PDF set MRST2004FF3 is used with 
m b = 4.3 GeV and (i r = (if = [1q = \JQ 2 + 4m 2 . The fraction of beauty quarks producing 
D*± mesons is taken as f(b -> D*) = 17.3 ± 2.0% fl3TU. 



4 Event Selection and Signal Extraction 

DIS events are selected by requiring a compact electromagnetic cluster in either the LAr or 
SpaCal calorimeters, which is taken to be the energy deposit of the scattered electron. The 
cluster has to be associated to a track reconstructed in the CTD. The events are triggered by 
either a coincidence of a SpaCal cluster and a signal from the CJC, or by the presence of 
a LAr cluster and a signal from the proportional chambers. The hadronic final state (HFS) 
particles are reconstructed using a combination of tracks and calorimeter deposits in an energy 
flow algorithm ||32l which avoids double-counting. The event kinematics including the photon 
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virtuality Q 2 , the Bjorken scaling variable x and the inelasticity variable y are reconstructed 
with the eS method [33], which uses the scattered electron and the HFS. The measurement is 
performed in the kinematic region 100 < Q 2 < 1000 GeV 2 and 0.02 < y < 0.7. 

The D* ± mesons from the decays D* ± (2010) D°(1865)tt^ ou , (K T 7r ± )vr^ are 
reconstructed using the tracks in the CTD. The branching ratio for this channel amounts to 
B = 2.63 ± 0.04% [|34l . The invariant mass of the Kn combination is required to satisfy 
\m(Kir) - m(D°)\ < 80MeV where m(D°) = 1864.84 MeV 01. The decay angle 6* of the 
kaon in the rest frame of the D° is restricted to cos 9* > —0.7, in order to reduce the back- 
ground, which strongly increases towards cos(#*) = — 1 as opposed to the D°, which decays 
isotropically. To further reduce the combinatorial background, a Q 2 -dependent cut on the D*^ 
transverse momentum, p T (D*) / GeV > (3 • [log(Q 2 / GeV 2 ) — 2] + 2), is applied. This criterion 
accounts for the increasing transverse momentum of the hadronic final state with rising Q 2 . 

The D*^ candidates in the pseudorapidity range |^(D*)| < 1.5 are selected using the mass 
difference method 11351 . In Fig. [[J a) the distribution of the mass difference Am = m{K^T\) — 
m(Kn) is shown for the selected data sample. A clear peak is observed around the nominal 
mass difference of 145.4 MeV. Wrong charge K ± n ± n T combinations with K ± tt ± pairs in the 
accepted D° mass range are used to describe the combinatorial background. 

The number of D*^ 1 mesons is determined in each analysis bin from a simultaneous fit to the 
signal and the background distributions. The Crystal Ball function [|36l is used for the signal 
description and the Granet parametrisation ll37l for the background. Several fit parameters in the 
single and double-differential distributions are fixed using the full data sample and the Monte 
Carlo predictions ll38~l . 

The cross section presented in this paper corresponds to the kinematic range summarised in 
Table |2l The p T (D*) and t](D*) range is chosen to be the same as in previous HI analyses HI 
at lower Q 2 . The Monte Carlo simulation is used for the extrapolation down to px(D*) = 
1.5 GeV. This extrapolation typically leads to a 15% increase in the cross section. With all the 
selection cuts, the average acceptance amounts to around 30%. 

Photon virtuality Q 2 100 < Q 2 < 1000 GeV 2 

Inelasticity y 0.02 < y < 0.7 

Pseudorapidity of D* ± -1.5 < rj{D*) < 1.5 

Transverse momentum of D*^ 1 p T (D*) > 1.5 GeV 

Table 2: Definition of the kinematic range of the present analysis. 



The inclusive D* 7 ^ production cross section is studied differentially in the kinematic vari- 
ables Q 2 , x, px(D*), i](D*) and the D*^ inelasticity z(D*), which corresponds to the fraction 
of the virtual photon momentum carried by the D* ± meson. The D* ± inelasticity is determined 
as z(D*) = P ■ pr,* I P ■ q = (E — p z )o* / ^yE e , where E e is the energy of the incoming electron 
and P, q and pr>* denote the four-momenta of the incoming proton, the exchanged photon and 
the D*^ 1 meson, respectively. The cross section for D*^ meson production is calculated from 
the observed number of D* ± candidates N D *±, according to: 

a vls (e + P - e+D*X) = ^ ' { )~"\ , (1) 

l~>int • -D • e • ( 1 + rad ) 
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where e is the reconstruction efficiency, r the contribution from reflections, C int the integrated 
luminosity, B the branching ratio and 5 ra d denotes the radiative corrections. 

The reconstruction efficiency accounts for the trigger efficiency and the detector acceptance 
and is determined using the Monte Carlo simulation. For this purpose charm DIS events are 
generated using both the RAPGAP and CASCADE programs and the average efficiency is used. 
For the efficiency determination, RAPGAP is reweighted in Q 2 and CASCADE is reweighted 
in pt(D*) in order to optimise the data description. The kinematic distributions of the D*^ 
candidates compared with the reweighted Monte Carlo predictions are shown in Fig.Qlb)-(d). 

The contribution r of reflections in the D° mass window from D° decay channels other than 
that considered in this analysis is estimated using the Monte Carlo simulation. This contribution 
amounts to r = (4.4 ± 0.5)% independently of the D*^ transverse momentum. The radiative 
corrections 5 ra d are determined using RAPGAP interfaced to HERACLES 4.1 ||39l and amount 
to 3% on average. The photoproduction background estimated using data If38l is not subtracted, 
but does not exceed 2.7%. The fraction of D*^ mesons originating from bb events is estimated 
as described in section [3j It amounts to 4% on average and is included by definition in the 
inclusive D*^ cross section. However, for the extraction of F 2 CC , the predicted contribution from 
bb production is subtracted from the data. 

5 Systematic Uncertainties 

The systematic uncertainties are estimated by varying the input parameters to the Monte Carlo 
simulations within the experimental precision at the reconstructed level or the range allowed by 
the theoretical models at the generator level. The following correlated uncertainties are taken 
into account: 

• The uncertainty on the hadronic energy scale is propagated to the measurement by chang- 
ing the hadronic energy by ±2% (±3%) for events where the scattered electron is detected 
in the LAr (SpaCal) calorimeter. The uncertainty due to the scattered electron measure- 
ment is estimated by varying the electron energy by ±1% and the polar angle by ±3 mrad, 
respectively. 

• The trigger efficiency, luminosity and D* — > Ktttc branching ratio are known with uncer- 
tainties of 1%, 3.2% and 1.5%, respectively. An uncertainty of 1.2% on the cross-section 
measurement arises due to the uncertainty on the photoproduction background. 

• The uncertainty on the reconstruction efficiency is taken as half of the difference between 
the two simulations, RAPGAP and CASCADE. This also covers the uncertainty on the 
extrapolation to pt(D*) = 1.5 GeV. The uncertainty in the efficiency determination due 
to the charm fragmentation model is estimated by varying the Kartvelishvili parameter 
a within its error as described in section [3] The uncertainty due to the choice of PDFs 
is estimated by using the CTEQ6L(LO) [|40]| parton densities in RAPGAP and the A2 
set EQ in CASCADE as alternatives. 

The following uncorrelated systematic uncertainties are accounted for: 
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• The signal shape and the invariant mass resolutions of the data are not fully reproduced 
by the Monte Carlo simulation. The errors on the signal extraction are determined 
by varying the fit parameters within their uncertainties. The fraction of events outside the 
D° mass window is determined using the Monte Carlo simulation. Half of this fraction is 
taken as a systematic error to account for the uncertainty on the D° mass resolution. 

• An uncertainty of 0.5% is assigned to the contribution from reflections to account for a 
possible p T dependence. The uncertainty of the QED radiative corrections is 1.5%. 

The following uncertainties are treated as partly correlated: The charged particle reconstruction 
uncertainty of 2.17%, which translates to 6.5% per D*^ and the uncertainty on the electron 
track-cluster matching of 2%. The above uncertainties are added in quadrature to derive the 
experimental systematic error. 

The theoretical uncertainties on the HVQDIS prediction are estimated by varying the in- 
put parameters as follows. The charm mass is varied from 1.3 to 1.6 GeV. The factorisation 
and renormalisation scales jjf = jj r are varied simultaneously from 0.5/i to 2/j . The frag- 
mentation parameter is varied within its error as described in section [3] The parton density set 
CTEQ5F3 ll42l is used as an alternative to MRST2004FF3. The resulting uncertainties, together 
with the error on f(c — > D*), are added in quadrature and are correlated between the bins. The 
uncertainties on the ZMFVNS prediction [12J are estimated by variation of the renormalisation 
and factorisation scales simultaneously from 0.5/io to 2/x . 



6 D* ± Production Cross Section 

The total inclusive cross section for D*^ 1 production in the phase space covered in this analysis 
(Table [2]) is measured to be: 

°vis(e + p -> e+D^X) = 225 ± 14(stat.) ± 27(syst.) pb . 

The corresponding predictions from RAPGAP, CASCADE and HVQDIS amount to 322 pb, 
279 pb, and 241^ pb, respectively, including the bb contribution. In Fig. [2] and Table [3] differ- 
ential cross sections are presented as a function of the DIS kinematic variables x and Q 2 and 
as a function of the D* variables pr(D*), r/(D*) and z(D*) . The data are compared to the ex- 
pectations from the HVQDIS calculation and from the RAPGAP and CASCADE Monte Carlo 
simulations. Neither Monte Carlo simulation describes the shape and normalisation of the D*^ 
kinematic distributions well, in contrast to the measurement [ 1 J at lower Q 2 . The HVQDIS 
calculation agrees with the data within the theoretical uncertainties. 

In Fig. [3] and Table |4] the double differential cross sections are shown as a function of y for 
different bins in Q 2 . The data are compared to the expectations of the HVQDIS calculation as 
well as to the RAPGAP and CASCADE simulations. HVQDIS describes the data well. Except 
for the first (Q 2 ,y) bin, the same holds for CASCADE. RAPGAP significantly overestimates 
the visible cross section. 
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The data are also compared to the ZMVFNS prediction [12J. This calculation has an in- 
trinsic limitation on the transverse momentum in the photon-proton center of mass frame, 
namely p^{D*) > 2 GeV. Therefore the same additional cut is applied to the data and the 
cross section is determined for the corresponding phase space. In Fig. |4]the D*^ cross sec- 
tions are shown as a function of pj,(D*), px(D*), rj(D*) and Q 2 , together with the ZMVFNS 
and HVQDIS calculations. The ZMVFNS prediction fails to describe the data, while HVQDIS 
agrees well with the data. 



7 Extraction of Ff 

The charm contribution -F 2 cc (x, Q 2 ) to the inclusive proton structure function F 2 is defined by 
the expression for the single photon exchange cross section for charm production: 



d 2 a c ~ c 2na 



cm 



<lnlQ2 Qij . ([l + (l-y) 2 ]Ff(x,Q 2 )-y 2 Ff(x,Q 2 )) , (2) 

where a em is the electromagnetic coupling constant. Weak interaction effects are neglected. 

The contribution from the structure function Fjf amounts to at most 3% lfi~3l in the present 
phase space and is neglected. The visible inclusive cross sections a^(y, Q 2 ) in bins of y 
and Q 2 are converted to a bin centre corrected F| c ((x), (Q 2 ) ) in the framework of a particular 
model using the relation: 



(Q 2 )) = a to< )Q n l ■ Ff theo ((x), (Q 2 )) , (3) 

[y,Q 



where and F| c theo are the theoretical predictions from the model under consideration. 
As in previous publications (TJ [2j the HVQDIS program and another program lfi~3l are used to 
calculate these quantities at NLO. CASCADE is not used for an F 2 CC extraction since it does not 
agree with the data (Fig. [2]). 

The model uncertainties on the measurement of F 2 CC are estimated by varying the HVQDIS 
parameters as described in section[5] The variations are made simultaneously in the calculation 
of the visible cross sections and in the prediction for F 2 CC . The total model uncertainties 
amount to 1 — 7% and are dominated by the variation of the renormalisation and factorisation 
scales. The central values of F| c with experimental and model uncertainties are summarised in 
Table [5] The fraction of the total cross section in the visible phase space, as predicted by 

a fy Q2\theo 

HVQDIS and given by < n2 iXh%o » is a ls° quoted and varies between 0.4 and 0.7. 

In Fig.[5]F 2 cc is shown as a function of x for different values of Q 2 . The F 2 CC values are con- 
sistent with those obtained in an inclusive track measurement using the HI vertex detector infor- 
mation 0. The expectation from the recent PDF fit to inclusive DIS data, HI PDF2009 ll44ll . 
tends to overestimate the data. In Fig. [5jb) the measurements are compared to the massive 
FFNS calculation at NLO O and NNLO @f0 and to the GMVFNS predictions at NLO and 
NNLO ll45l l46ll . The FFNS predictions agree well with the data over the full kinematic re- 
gion investigated. The expectations for F 2 CC from a global fit in the GMVFNS at NLO tend to 
overestimate the data. At NNLO the GMVFNS prediction agrees better with the data. 
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8 Conclusions 



The cross section for D*^ meson production is measured in the phase space 100 < Q 2 < 
1000 GeV 2 and 0.02 < y < 0.7. Single and double differential cross sections are compared 
to Monte Carlo simulations and the predictions of NLO calculations in massive and massless 
schemes. The data have a typical precision of 20%. 

In the measured domain the RAPGAP and CASCADE simulations do not provide good 
a description of the D*^ kinematics. The double-differential cross section d 2 a/dy dQ 2 is de- 
scribed well by CASCADE, while RAPGAP overestimates the cross section at high Q 2 . The 
NLO FENS calculation HVQDIS agrees with the data well, while the calculation based on 
ZMVFNS fails to describe the data. 

The charm contribution F| c to the proton structure function F 2 is determined. HVQDIS 
is used for extrapolation of the visible D*^ cross sections to the full phase space in p T (D*) 
and T](D*). The model uncertainties are found to be small in the kinematic region studied. The 
data are compared to QCD predictions at NLO in the FFNS scheme and to the CASCADE 
implementation of the CCFM model as well as to the expectations from global fit analyses, 
using GMFVNS implementations at NLO and NNLO. Both FFNS and CASCADE describe the 
measurement well. The data indicate that the NLO FFNS provides the best description of D* 
production and of F 2 CC in the kinematic region of the analysis. 
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Table 3: Single differential cross sections for D* ± production in bins of Q 2 , x and the meson 
kinematics, pr(D*), rj(D*) and z(D*), as measured in the visible range defined in Tabled The 
central values of the cross section are listed together with relative statistical (S sta t), uncorrelated 
(duncor) and correlated (5 cor ) systematic uncertainties. 
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Table 4: Double differential cross sections for production in bins of Q 2 and y as measured 
in the visible range defined in Tabled The central values of the cross section are listed together 
with relative statistical (S stat ), uncorrelated (5 unc ) and correlated (S cor ) systematic uncertainties. 
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Table 5: The measured values and relative errors for the charm contribution to the proton struc- 
ture function F 2 CC . Relative statistical, correlated and uncorrelated experimental systematic as 
well as model uncertainties are listed. The fractions of the total D* ± cross section in the visible 
phase space as predicted by HVQDIS are also given. 
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Figure 1: a) Distribution of Am = m^Kim) — m(Kn) for D*^ 1 candidates (K^ir^irf) and 
for wrong charge combinations {K ± ti ± -kJ) in the accepted D° mass window. The fit func- 
tion is also shown. Comparisons at the detector level between the D*^ data sample and the 
reweighted Monte Carlo models are presented. Background-subtracted distributions are shown 
as a function of Q 2 (b), p T (D*) (c) and rj(D*) (d). 
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Figure 2: Differential cross sections for inclusive D*^ 1 meson production as a function of 
Pt(D*), i](D*), z(D*), Q 2 and x. The inner error bars indicate the statistical uncertainties, 
the outer error bars show the statistical and systematic uncertainties added in quadrature. The 
expectations of CASCADE (dashed line) and RAPGAP (solid line) are obtained using the pa- 
rameters as described in section [3] The band of the HVQDIS prediction (shaded) is obtained 
using the parameter variation described in section[5] The ratio R = cr theory /& data is also shown. 
In the case of HVQDIS the theoretical uncertainties are taken into account. The inner error bars 
on the data points at R = 1 display the relative statistical errors, and the outer error bars show 
the relative statistical and systematic uncertainties added in quadrature. 
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Figure 3: Double-differential cross sections for production as a function of y in different 
Q 2 bins. For the purpose of presentation the cross sections are multiplied by Q 4 . The data 
(closed symbols) are shown with the statistical (inner error bars) and total (full error bars) 
uncertainties. Predictions from the RAPGAP (solid line) and CASCADE (dashed line) Monte 
Carlo simulations and the HVQDIS NLO calculation (shaded area) are also shown. 
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Figure 4: Differential cross sections for inclusive D*^ 1 meson production as a function of 
p T (D*), r](D*), p* T (D*), Q 2 and x as measured for p* T (D*) > 2 GeV. The inner error bars 
indicate the statistical uncertainties, the outer error bars show the statistical and systematic un- 
certainties added in quadrature. The expectation of HVQDIS (shaded band) is obtained using 
the parameter variation described in section [5l The prediction in ZMVFNS is represented by 
the hatched band where the uncertainty originates from the scale variation. 
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Figure 5: The charm contribution F| c to the proton structure function. The data (closed sym- 
bols) are shown with statistical (inner error bars) and total (full error bars) uncertainties. In a) 
the data are compared to the HI measurement of F 2 CC using secondary vertex information (open 
symbols) 0, where measurements at Q 2 = 300 GeV 2 are shifted to Q 2 = 400 GeV 2 using the 
NLO calculation [13J. The result of the PDF fit H1PDF2009 (shaded band) is also shown. The 
uncertainty band accounts for experimental, model and parametrisation uncertainties Il44l . In 
b) the data are compared to the QCD predictions from the NLO calculation [fT3l in FFNS (light 
thick solid line). The predictions from the global PDF fits MSTW08 at NLO (dashed) and 
NNLO (dark solid) as well as the results of the ABKM fit gH at NNLO in FFNS (dotted) and 
GMVFNS (dashed-dotted) are also shown. 
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